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ABSTRACT 

The effects of relativistic expansion on the late-time supernova light curves 
are investigated analytically, and a correction term to the (quasi-)exponential 
decay is obtained by expanding the observed flux in terms of j3, where f3 is the 
maximum velocity of the ejecta divided by the speed of light c. It is shown 
that the Doppler effect brightens the light curve owing to the delayed decay of 
radioactive nuclei as well as to the Lorentz boosting of the photon energies. The 
leading correction term is quadratic in (3, thus being proportional to E k /(M e jC 2 ), 
where and M e j are the kinetic energy of explosion and the ejecta mass. It 
is also shown that the correction term evolves as a quadratic function of time 
since the explosion. The relativistic effect is negligibly small at early phases, but 
becomes of considerable size at late phases. In particular, for supernove having a 
very large energy(hypernova) or exploding in a jet-like or whatever non-spherical 
geometry, 56 Ni is likely to be boosted to higher velocities and then we might see 
an appreciable change in flux. However, the actual size of deviation from the 
(quasi-)exponential decay will be uncertain, depending on other possible effects 
such as ionization freeze-out and contributions from other energy sources that 
power the light curve. 

Subject headings: radiation transfer - sup ernovae: general 



1. Introduction 

It has long been recognized that the radiative transfer in supernova(SN) ejecta should 
be treated relativistically to account for the high velocities achieved in its outermost layers. 
However, it is only in the last decade that the first numerical codes for radiative transfer 
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that take relativistic effects into account were developed and applied to transfer problems in 
SN ejecta. At early times of explosion, the energy of photons emitted from the photosphere 
of SN ejecta is greatly enhanced by the Doppler effect. As time goes on and the photosphere 
recedes towards the center of ejecta, the degree of the enhancement decreases, and eventually 
the resultant light curve(LC) is expected to follow an (quasi-)exponential deposition curve 
due to radioactive decays. However, there are several factors that cause the late-time LCs 
deviate from the radioactive decay curve. One possibility is the flattening of LCs owing to 
the so-called ionization freeze-out effect, as first pointed out by Fransson & Kozma (1993) 
for the late-time LC of SN 1987A. Another one is the signatures of contributions from other 
energy sources such as possible pulsar activity and circumstellar interaction as have been 
discussed by previous works. 

In this Letter, I would like to point out that late-time SN LCs may show a deviation 
from radioactive decay curves as a result of the relativistic Doppler effect. In particular, 
for very energetic SNe or what are called hypernovae(HNe), the relativistic effect would be 
very important since the maximum velocity near the surface of ejecta reaches a significant 
fraction of the speed of light. It will be shown that the Doppler effect makes a LC brighter 
due to the Lorentz boosting, but at the same time the light-traveling-time effect and a pure- 
relativistic effect delay the decay of radioactive nuclei, thus tending to make the LC even 
brighter at late phases. The net result is determined by the sum of these two effects. Late- 
time LCs provide direct information to determine the 56 Ni masses ejected by SNe. Therefore, 
its precise determination is of great importance for studies of SN explosion mechanism and 
chemical evolution of galaxies. In this Letter, we study the relativistic effect on the late-time 
LCs and estimate the size of the correction to the (quasi-)exponential decay for ordinary SNe 
and HNe. For simplicity, we use the approximation that energy input from radioactivity is 
emitted as optical photons on the spot and carried away from the SN ejecta free of absorption. 



2. Light Curves 

We consider a spherically symmetric SN ejecta that expands homologously, so that its 
radius is given by R(t) = Rq + f3ct as a function of time t after explosion, where Ro is 
the initial radius, (3c is the surface velocity, and c is the speed of light. Let us introduce 
the radial and azimuthal coordinates, x and fi, as shown in Figure 1, to designate volume 
elements of ejecta, such that x runs from to 1 outwards and that /x is the cosine of angle 
9 made between radial direction and the line of sight. 

The flux of photons received by an observer at time t is the sum of the emissions from 
the volume element at a retarded time t rct , dV = 2n R(t ret ) 3 x 2 dxdfj J , where t Tet is defined by 
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c(t — t TCt ) = d — R(t rct )x cos 9 

~ d - c(3fixt ret (1) 

and d is the distance between the SN and the observer. We used the fact that the initial 
hydrodynamical time scale is short compared to the elapsed time considered, Rq/(5c « t. 
Solving equation (1) in terms of t rct , we obtain 



_ t-d/c 

tret " \^~x (2) 

Let us define the frequency-integrated emissivity j as the total energy of photons emitted 
per unit volume, unit time, and unit solid angle. Then, the emissivity in the rest frame in 
the direction toward the observer j(x, /i, t ret ) is related to that of the comoving frame j c such 
that 



J {X, ft, tret) = -gTj T. 7g , 3) 

71(1 -(3fix) 3 

where j x = (1 — p 2 x 2 )^ 1 ^ 2 and t rc t, c = trct/lx is the retarded time in the comoving frame 
(e.g.,Rybicki & Lightman 1979). We assume that the radiative loss from the ejecta is bal- 
anced by the energy input due to radioactive decay and that the emissivity in the comoving 
frame is isotropic. Then, j c is written as 



j c (x,/i,t rctjC ) = — /e , (4) 

47T r 

where n c is the number density of radioactive nuclei in the comoving frame, r is its decay 
time, e is the energy available per decay, and / is the deposition fraction. The total number 
of radioactive nuclei in a volume element AV^, r y x n c AV, obeys the decay law given by 



or 



(5) 
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dt ret 'JxT 

which has a solution 



" n c AV) = -^, (6) 



^rct Q ^rct 



n c (x,t rct ) = ^77^(*,0)e 7,r = ) n c{X}0 )e 7,r 

2 ^ret 

i?o \ Ji(x, 0) 



^R(t ret )J 7a; 

Using equation (7) in equation (4), we obtain 



e 7*t. ( 7 ) 



tr 



3 * 

i?o \ fen(x,0) 



In equating the energy emitted in a time <it ret by a volume element at x, fi and time 
t re t, dV(x, t ret ), with the corresponding energy passing through a normal area AS at the 
observer during time dt, ASdF(t)dt, we have 



dF(t)dtAS = j(x, n, t Tet )dt Tet dV(x, fx, tret) An, (9) 

where dF(t) is the differential flux corresponding to the volume element and AQ is the solid 
angle subtended for the area AS by the volume element, and thus AS* = d 2 AQ. Adding up 
contributions by all the volume elements, we have the observed flux at time t, 



F(t) = j 2 J^j(x,^t ret )dV 

1 fN tot e f 1 2 1 f 1 1 

/ dxSx n(x) x - / da—- — rrexp 



And 2 r 



(t - d/c) It 
7 X (1 - pLix) 



.(10) 



Here, n(x) = n(x, 0)/n is the dimensionless number density of the radioactive nuclei 
and iVtot = (47r/3)-RQn is its total number at t — 0, where n is the average number density 
defined as 
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/ n(x, 0)x 2 dx x 
n = ^ j = 3 / n(x, 0)x 2 dx, (11) 



/ 

Jo 



x 2 dx 

and thus n is normalized such that 



f 3x 2 h(x)dx = 1. (12) 
Jo 

In deriving equation (10), it is assumed that / is constant throughout the ejecta for simplicity. 
This holds for 56 Co decay at sufficiently late times when the deposition is primarily due to 
e + decay, and even for general cases if only the effect of a varying / is incorporated into n. 
Using the variable y = 7~ 1 (1 — fifix)' 1 instead of //, the equation (10) is rewritten as 



F(t) = 47r ^ 2 fN ™ e f dx3x 2 h(x)m x, X), (13) 



where 



1 



Zfxpx J a 



b(P,x) 



f(J3,x,X) = —— y 2 e- Xy dy, (14) 



09,*) 



X = — ft ) , and a((3, x), b((3, x) are given by 



T \ C 



a \P' x ) = h I a ,b((3i x )=[-, 7T\ ' ( 15 ) 



respectively. Note that b(/3,x) = a(—/3,x) so that f(/3,x,X) is an even function of (5. 
The integration in equation (14) can be readily done to yield 



m *' X) = 2^fx ' (16) 
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with 



g(y) = J y 2 e- Xy d y = ~e~ x y { X V + 2X v + 2 ) • ( 17 ) 

Expanding g(a(/3,x)) and g(b(/3,x)) in equation (16) in terms of f3, we obtain 



f((3,x,X) = g'(l) + 



"(1) , ^ (3) (1) 



+ 



6 



(3 2 x 2 + 0(/T) = e 



-x 



1 + 



X 2 - 7X + 
6 



(3 2 x 2 + 0(/3 4 ) 



.(18) 



Using equation (18) in equation (13) and retaining terms up to the second order in (3, 
we get 



Fit) - 1 /^tote x 



1 + 



X 2 - 7X + 



(19) 



where < > is the average of square velocity of radioactive nuclei given by 



< p 2 x >= [ dx3x 2 h(x)p 2 x 2 . 
Jo 



(20) 



3. Size of the Relativistic Correction 

The leading term in equation (19) corresponds to an (quasi-)exponential decay of the 
flux that is expected for the exponential decay of radioactive energy sources. It is quasi- 
exponential, because, in general, / changes with time. The next term is a relativistic correc- 
tion, which is of second order in f3 and is generally very small. The correction is quadratic 
in t' = t - d/c. It takes a value (4/3) < f3 2 > at t' — and has a minimum -(17/24) < 0% > 
at t' = 3.5t. At late times when t' >> r, the correction term becomes quite large. 

To calculate < (3 2 >, we need to know the density structure of the ejecta. Here we 
assume a density profile 



p{x) 



< x < x c , 

X c < X < 1, 



(21) 
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which is a good approximation to the actual SN ejecta and greatly simplifies our analysis. 
This profile is characterized by the core radius x c and the index of the power-law p. For 
given x c and p, p c is determined to be 



P - 3 3M ej 
c px\ — 3xc 4:%R 3 

In addition, the maximum velocity achieved in SN ejecta, /3, is given by 



of IQgk p-5 yEg-3xg y /2 

P V 3M ejC 2 p-3 px*-5a%J ' 1 J 

where £"k and M e j are the kinetic energy of explosion and the ejecta mass, respectively. Using 
the fact that x c ~ 0.1 and p is as large as ~ 8 — 10 for typical cases, the equation (22) reduces 
to 



which gives (3 ~ O.33(E k /lO 51 erg) 1 / 2 (M cj /M )- 1 / 2 (a; c /O.l)- 1 for p = 8. 

If radioactive nucleus 56 Ni is homogeneously mixed within the ejecta and thus n(x) has 
the same profile as p(x), we find from equation (20) 



< ^ >= -.P / 2 4) 

5 p — 5 pxj? — 3xc 

From equations (22) and (24), we obtain < 01 >= 2_E k /(M cj c 2 ). Then, the correction 
term in the bracket of equation (19) turns out to be 



*Iz"LH> < A > = ni x 10 -3£!^ (^ (^ ) . (25 ) 



u 6 V10 51 erg/ V^o, 

It is seen from equation (25) that the size of the relativistic correction depends only on 
M e j and E k , irrespective of the density profile of the ejecta in the homogeneous-mixing case. 
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We use the fact that radioactivity is dominantly mediated by 56 Co decay of the decay 
chain 56 Ni^ 56 Co^ 56 Fe, and thus r = r Co - r Ni = 111.3 - 8.8 = 102.5 days and e = e Co , the 
energy available per decay of 56 Co, effectively. Then, it is found that the correction term is 
only 0.16, 0.88, and 3.7 percent of the leading term even at t' =2, 3, and 5 years, respectively, 
for the case of Type la SN(SN la), which has a relatively large expansion velocity because 
of its small mass M C j ~ M & and canonical explosion energy E k = 10 51 erg. For Type II 
SNe(SNe II), the ejecta mass is much larger, e.g., M ej ~ 1OM , so that the correction will 
be much smaller. In the actual SN ejecta, 56 Ni is more centrally concentrated than the 
homogeneous-mixing case now considered. Thus, the above numbers are likely to be an 
overestimate. Therefore, it is concluded that, for ordinary SNe, the relativistic effect would 
be very small compared to other possible effects such as ionization freeze-out, circumstellar 
interaction, and contributions from radioactivity of other than 56 Co. 



4. The Case of SN 1998bw 

HNe have larger explosion energies than ordinary SNe by an order of magnitude. Thus, 
we expect that the relativistic effect will be more important for HNe. As an example of HNe, 
we take a model for a hypernova SN 1998bw, C0138 (Iwamoto et al. 1998; Nakamura et al. 
2001), which has M cj ~ 1OM and E k ~ (20 - 50) x 10 51 erg. We found that the density 
distribution of C0138 is well reproduced by a profile given by equation (21) with p = 10 and 
x c = 0.1. Figure 2 compares the analytic fit (dotted line) and the result of hydrodynamical 
calculation (solid line) for n of C0138 with E k = 30 x 10 51 erg (Iwamoto et al. 1998) and 
with homogeneous 56 Ni distribution. The value of p is in good agreement with results by 
Matzner & McKee (1999) obtained for cases of compact SN progenitors. 

Then, from equation (25), we find that the change of the flux will be still relatively small 
at early times, but it grows considerably at late times. It is only 1.8, 2.6, and 4.4 percent 
at t' — 3 years, but becomes 7.4, 11.1, and 18.5 percent at t' — 5 years for cases of E k = 
20, 30, and 50 x 10 51 erg, respectively. If the explosion occurs in a jet-like (Nagataki 2000) 
or whatever non-spherical geometry, 56 Ni is likely to be boosted to higher velocities and 
then we might see an appreciable change in flux. In late-time spectra of SN 1998bw, it was 
observed that iron has a larger velocity than oxygen, indicating that such a non-spherical 
56 Ni distribution may be realized in the ejecta of SN 1998bw (Mazzali et al. 2001). 

Even for spherically symmetric models, it is reported that the ejecta of HNe may become 
trans-relativistic with its outer envelope having a mean velocity 7 = (1 — f3 2 )~ l l 2 ~ 2 {(3 ~ 
a/3/2) and a flatter density gradient of p ~ 4 (Tan, Matzner, & McKee 2001). Therefore, 
to make a detailed comparison between the observed late-time LCs and those predicted by 
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non-spherical HN models, it is recommended that we calculate model LCs by taking the 
relativistic effect into account. Note that for large values of (3 the expansion series in terms 
of f3 (equation 18) may not converge well. However, the analytic solution with equations 
(16) and (17) is exact, and thus we can use it to study such cases with large f3. 
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Fig. 1. — Schematic picture showing the position of a volume element of SN ejecta at a 
retarded time t ret . 
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Fig. 2. — Normalized number density of 56 Ni, n as a function of x, in the case of homogeneous 
mixing for model C0138 of SN 1998bw. The solid line is the result of hydrodynamical 
calculation (Iwamoto et al. 1998), while the dashed line is an analytic fit using equation (21) 
with p = 10 and x c — 0.1. 



